Persistence of a Rouse polymer chain under transverse shear flow 
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We consider a single Rouse polymer chain in two dimensions in presence of a transverse shear 
flow along the x direction and calculate the persistence probability Po(i) that the x coordinate of 
a bead in the bulk of the chain does not return to its initial position up to time t. We show that 
the persistence decays at late times as a power law, Po{t) ~ t~^ with a nontrivial exponent 6. 
The analytical estimate of S = 0.359... obtained using an independent interval approximation is in 
excellent agreement with the numerical value 9 « 0.360 ± 0.001. 
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I. INTRODUCTION 



Polymer dynamics plays a central role in material sci- 
ence and biology. In particular, dynamics of an indi- 
vidual flexible or semi-flexible polymer under a suitable 
shear force has been of great interest P, [3, H, 0, H, d, 0| ■ 
Shear force comes into play when a fluid flows past a 
surface. Substantial effort has been undertaken to in- 
vestigate the motion of polymers in shear field. Pre- 
viously, studies were done on bulk samples using light- 
scattering and birefringence experiments. Recently, the 
dynamics of a single polymer has also been investigated 
using video-microscopy. Under a shear stress, such a 
polymer shows tumbling in addition to a longitudinal 
stretching [1, If one takes a tethered polymer, whose 
one point is made immobile, then tumbling leads to a 
cyclic motion of the spatially constrained polymer about 
a mean position |^]. Also statistics of polymer orien- 
tation angles has been of interest [1, 0] ■ These proper- 
ties studied both experimentally and theoretically can be 
classified as long-time transport phenomena. In contrast, 
in this paper we explore the persistence or the survival 
probability behaviour of a flexible polymer chain under 
transverse shear flow within the paradigm of the simple 
Rouse model where the polymer chain consists of beads 
or monomers connected by harmonic springs 0] . We will 
show that even in this simple model, the persistence at 
late times decays as a power law characterized by a non- 
trivial exponent. 

The survival/persistence probability Po(i) that a 
stochastic process X{t) does not cross zero upto time 
t is a quantity of long standing interest in probability 
theory and with many practical applications @. The 
derivative F{t) — —dPo{t)/dt is the first-passage proba- 
bility [lOj. In many nonequilibrium many body systems, 
the persistence has been found to decay as a power law 
at late times, Po{t) ^ t^^ . The exponent 6 is called the 
persistence exponent and has been a subject of much 
theoretical, numerical and experimental studies in re- 
cent times [llj |. The exponent 9 is often nontrivial and 
is generally hard to calculate analytically even in sim- 



ple systems such as the linear diffusion equation start- 
ing from random initial conditions The reason for 
this difficulty can be traced back to the fact that the 
spatial interactions in these extended systems makes the 
local stochastic field X{t) a 'non-Markovian' process in 
time 

In this paper we study the persistence properties of a 
Rouse chain in 2-dimensions in presence of a transverse 
shear velocity field which is non-random. We show that 
the persistence probability in this system decays at late 
times as a power law with a nontrivial persistence expo- 
nent 9 ~ 0.36 that we compute numerically as well as an- 
alytically within an independent interval approximation 
(IIA). We note, that the current problem is in contrast to 
similar problems considered in "random" flow fields ear- 
lier. For example, for a Rouse chain ^ of infinite length, 
the transport properties [H, [3, [lB| and the persistence 
properties jl6| in a quenched random velocity flow field 
have been studied. 

The paper is organized as follows. In Section II, we de- 
fine the model precisely and summarize our main results. 
In Section-Ill we present exact calculations of the two- 
time correlation functions in our model. These results 
are used next in Section IV to calculate the persistence 
exponent analytically within the IIA. Sections V and VI 
describe details of the numerical methods and finally we 
conclude in Section VII. 



II. THE MODEL AND MAIN RESULTS 

We consider a Rouse polymer chain embedded in a 
2-dimensional plane. The chain consists of beads con- 
nected by harmonic springs @. In addition, the chain 
is advected by shear velocity fiow field. Let [a;„(t), y„(t)] 
denote the coordinates of the n-th bead at time t which 
evolve with time according to the following equations of 
motion 



dyn -p , 

—J- = r (a;„+i -I- a;„_i - 
at 



2y^)+m{n,t) (1) 
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where T denotes the strength of the harmonic interaction 
between nearest neighbour beads, r/i{n,t) and 772(71,, t) 
represent the thermal white noises along the y and x 
directions respectively that are uncorrelated. The trans- 
verse shear velocity field v{y) is linear 



v{y) = y. 



(3) 



For a finite chain with N beads, Eqs. ([T]) and ([2]) are 
valid only for the (iV — 2) interior beads. The two bound- 
ary beads will have slightly different equations of motion. 
However, for an infinitely large chain {N — > 00), the 
translational invariance along the length of the chain is 
restored since the boundary conditions become irrelevant 
for late time dynamics. Since we are mostly interested in 
the late time properties, one can make further simplifi- 
cations by replacing the discrete index n of the beads by 
a continuous variable s and subsequently replace the dis- 
crete Laplacian by a continuous second derivative along 
the s direction. The coarse grained versions of the evo- 
lution equations Q then become 



dy{s,t) ^ dMs,t) , 

= r — — h?7i(s,i). 



dt 
dx{s, t) 
dt 



ds^ 



+ 2/(s,t). 



(4) 
(5) 



Note that we have also dropped the 7/2 (s,t) term in the 
second equation. This is simply because one can easily 
show that the noise term 772 (s,t) becomes insignificant 
compared to the shear force term y{s,t) at late times. 
Hence for late time asymptotic properties we can ignore 
the noise 772(5, t). 

In the absence of harmonic interactions (F = 0), the 
beads become independent and the coordinates of any 
(say the n-th) bead represents a two-dimensional Brown- 
ian walker in a shear flowfl7'|. Equivalently , in this limit, 
the X coordinate of the walker evolves as (Pxn/dt^ ~ 
r]i{n,t), i.e., it represents a randomly accelerated parti- 
cle. The persistence probability of the x-coordinate, i.e., 
the probability that the x coordinate does not cross zero 
up to time t is known to decay as ^ t~^l'^ [l^. Re- 
cently, the persistence of a single random walker for var- 
ious other deterministic velocity functions v(y) has also 
been studied Interestingly it has been shown 

that for all odd functions v(%j) survival probability de- 
cays as t^^l^ [19]. It turns out that the same t~^l'^ decay 
also holds in the case when w(y) is not a deterministic 
function, but represents a quenched random transverse 
velocity field with short-range correlations [2l|, [22] ■ This 
model of a single random walker in presence of a random 
transverse velocity field is known as the Matheron-de- 
Marsily modelfSs^ whose transport properties had been 
studied earlier extensively [13], but the studies of persis- 
tence properties are relatively new[2l[. [2^. [25j. 

In this paper we study the persistence probability of 
the X coordinate of the 77,-th bead in the presence of har- 
monic interaction F 7^ 0. Due to the translational invari- 
ance along the length of the chain in the bulk, the per- 
sistence probability is independent of the label n of the 



bead. We also absorb the factor F by properly rescaling 
the time. Note that the continuum equation ([4]) for the 
y coordinate is precisely the Edwards- Wilkinson equa- 
tion of one dimensional interface [1^ and its persistence 
properties are known, both theoretically [l^, [1^ and also 
experimentally [2^. Here we focus on the x coordinate 
and define the persistence as follows 



Po(i) = Prob]a;(s,t') 7^ a;(s,0) 
for all t! : < t' < t], 



(6) 



i.e., Po{t) is the probability that the x coordinate of any 
bead does not return to its initial position within the 
time interval [0,t]. 

The initial conditions for the chain coordinates do not 
play any role in the persistence probability. This is due 
to the fact that the evolution equations are linear, so we 
can redefine the change in positions x{s, t) — x{s, 0) and 
y{s,t) — y{s,0) as the relevant coordinates which satisfy 
the same evolutions equations. Hence, for the evolution 
equations ^ and (O we can set the initial conditions 
x(s, 0) = and y{s,0) — without any loss of generali- 
ties. 

Our main results can be summarized as follows. We 
show that Po{t) ~ t~^ at late times t where the persis- 
tence exponent 9 has a nontrivial value. The numerical 
value 6 ~ 0.360 ± 0.001, is in excellent agreement with 
the analytical value 9 — 0.359... obtained within the IIA 
method. Thus as one switches on the harmonic interac- 
tion T ^ between the beads, the exponent 6 ~ 0.36 
increases from its value 9 = 1/4 for F = 0. Thus the 
X coordinate of a bead survives less in presence of har- 
monic interactions, i.e., the interaction enhances the re- 
turn probability. 



III. CALCULATION OF EXACT TWO-TIME 
CORRELATION FUNCTIONS 

The stochastic processes x(s, t) and y(s, t) evolving via 
Eqs. ((4|) and ((5|) are both Gaussian at late times since 
the evolution equations are linear. A Gaussian process is 
completely specified by its two-time correlation function. 
More detailed quantities such as the persistence prob- 
ability, in principle, is a complicated functional of the 
two-time correlation function. In this section we com- 
pute the two-time correlation functions exactly and use 
these functions later for computing the persistence prob- 
ability in Section IV. 

To begin with, we Fourier transform Eq. SI We de- 
fine y{k,t) — y{s,t) exp{—isk)ds, and f]i{k,t) — 
rii{s,t) exp{—isk)ds. This implies, 



dyik,t) 
dt 



-ey{k,t)+^\ik,t) 



(7) 



Assuming fiat initial condition (i.e., y{s,t) — 0), Eq. ([7]) 
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gives 



= exp(-fc2i) / rji{k,t')ex.p{k'^t')dt', (8) 
Jo 

which in turn impHes the correlation function 



(9) 



For Eq. defining x(k,t) — 

x{s,t) exp{—isk)ds and again considering flat 
initial condition (i.e., a;(s,0) — 0), we get 

S:{k,t)^eyiY>{-kH) [ y{k,t') ex.j>{kh')dt' , (10) 
Jo 

which further implies, 

{x{ki,ti)x{k2,t2)) = exp(-fci^ti - /C2^i2) x 

' dt' f dt"{y{ki,t')y{k2,t")) exp{kft' + kit"). 
Jo 

(11) 

Substituting Eq. ^ in Eq. HI]), we get 



where B is an unimportant constant. 

Note that due to the translational invariance in the 
bulk, the correlator of the process x{s, t) does not depend 
on the location s of the bead along the chain. Thus, for 
simplicity of notations, we can now drop the label s and 
consider x{t) as the relevant Gaussian process with the 
correlator C{ti,t2) = {x{ti)x{t2)) as given in Eq. pi)) . 
Clearly the process x{t) is non-stationary since its two- 
time correlator in Eq. (fT4|) depends on both ti and t2 
and not just on their difference. One can however define 
a logarithmic time T = In t and consider the normalized 

The survival or 



31| 



process X{T) = x{t) / {x{ty) in T 
no zero crossing probability is clearly the same for both 
the normalized process X{T) and the original unnormal- 
ized processs x(t). It then follows from Eq. p4)) that 
the autocorrelation function of this normalized Gaussian 
process X{T) is stationary in the T variable and is given 

by 



A{T) 



^Cih,h)C{t2,t2)' 



-^exp(^T) [exp(-r)Vl + exp(-r) 



-i(l + exp(-T)f/^ 



(l-cxp(-T)) 



5/2i 



(15) 



{x{ki,ti)x{k2,t2)) 

S{ki +k2) I dt' 



exp(- 







This form of the stationary autocorrelator will be used in 
the next section to compute the persistence probability. 



dt" 



t2-2 mmjt', t"))) - exp(-fci'(^i + tz)) jv. 



2ki' 



CALCULATION OF THE PERSISTENCE 
EXPONENT e 



(12) 



By inverting the Fourier transform above, we obtain the 
correlation function 



C{ti., 
ti 



t2) 



dt' 











dko 



exp(i(fci + k2)s) 

ti rt2 

dt / dt 
Jo 

1 



exp(- 



X{s,tl)x{s,t2)) 

dt" I dki 

oo J ~oo 

x{ki,t )x{k2,t )) 
dki 

-oo 

+ t2)) 1 -exp(-fci2(ti +<2 -?)) 



ki\t 
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We have thus mapped our problem to a Gaussian sta- 
tionary process in T = Int variable with a prescribed 
correlator A(T) and we want to calculate the probability 
Po{T) that the process does not cross zero up to time T. 
For a general correlator A{T), the computation of Po{T) 
is very hard [ll|, H^l- However, some general results 
are known for the late time behavior of Po{T). For ex- 
ample, it is known [ll|, [s^l that when A(T) decays 
faster than 1/T for large T, the persistence probability 
Po{T) decays exponentially, Po{T) ~ exp(— 6'T). Since, 
in our case, A{T) is Eq. p5|) decays faster than 1/T for 
large T, we expect Po{T) ^ exp{—9T). In terms of the 



,original time variable, t 



this would signify a power 



where t 
leads to. 



2 min(t',i"). After some algebra, Eq. ([T 



C{tut2 



B 



tlt2\/tl + t2 



{{ti+t2f/^ + \t2~t,\'>r'} 



5/21 



(14) 



law decay of the persistence Po{t) t~ for large t. Thus 
(13);he inverse decay rate 9 in the T variable is precisely the 
exponent of the algebraic decay in the real time t. 

While we were not able to compute the exponent 9 
exactly, one can obtain a very accurate analytical es- 
timate of 9 using the IIA method that was first used 
in the context of persistence in diffusion equation p^ . 
This method works reasonably well only for smooth Gaus- 
sian stationary processes. A process is smooth if A{T) = 
1 — aT^ + ■ ■ ■ for small T. In that case, the process has 
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a finite mean density p = ^J-A"{Q)Itt [335 of zero cross- 
ings. For our process, the correlator in Eq. (fTS]) can be 
expanded for small T 

A(r^O) = l-— t2 + ^t5/2, (16) 
16 V2 

indicating a = 15/16 and thus proving that the process 
is smooth. 

In the IIA, applicable only to smooth processes, one as- 
sumes that the intervals between successive zero crossings 
of a Gaussian stationary process are statistically indepen- 
dent. Within this approximation, one can then express 
the distribution P{T) of the intervals between successive 
zero crossings in terms of the correlation function A{T) 
in the Laplace space [l2| 



P{s) 



i-((r)/2),s[i-,si(.)] 



(17) 



l + {{T)l2)s[l^sA{s)] 

Here P{s) and A{s) are the Laplace transforms of P{T) 
and A{T) respectively and (T) = 1/p = 7r/^-A"(0) is 
the mean interval size. 

Using the exact expression of A[T) from Eq. (fTS)) . 
we can find P(T) from the above formula. The persis- 
tence probability Po{T) is simply related to the inter- 
val distribution [H, (PPo{T)/dT'^ = P{T)/{T). Since 
we expect Pq(T) to decay exponentially at late times T, 
i.e. Po{T) ^ exp(— 6'T), it follow that the interval dis- 
tribution P{T) will also have the same late time decay, 
P{T) ~ exp(— 6'T) with identical exponent 9. This means 
that the Laplace transform P{s) must have a simple pole 
at s = —9. In other words the denominator in Eq. pT|) 
must have a root at s = —6. Substituting s = —9, the 
denominator reads 

G{9) ^l-{e/2)i-'' 



2a 



1 + 9 A{-t 



(18) 



where we have put (T) = 7r/V2a with a = 15/16 and 9 
is given by the smallest positive root of G{9) = 0. 

To determine the root of G{9) — accurately, it is 
convenient to switch variables and define x = exp{—T), 
such that 

, \ 2x-l 

A{x) = 



1 



— a; -|- \/l — .T^ 



(19) 



and then Eq. (fTSl) becomes. 



2x~l 



l + \/l~ 



Solving Eq. ((20|) numerically gives 
6I11A = 0.359.... 



(20) 



(21) 



Thus the persistence probability Po{t) ~ decays al- 
gebraically for large time t with a nontrivial exponent, 
whose analytical value within the IIA is ^ha — 0.359... 
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FIG. 1: Persistence probability Po{t) versus t, simulated for 
the polymer chain length iV = 1000, 2000 and 10000. The 



fitted power law t 
the data. 



is shown by plotting a thick line above 



V. SIMULATION OF DISCRETISED 
LANGEVIN EQUATIONS 

In this section we describe simulation of the Rouse 
chain evolving via Eqs. ([T]) and ^ and further discre- 
tised in time t as 

Vnitm+l) = Vnitva) + At[j/„+i(tm) -|- y„-l(tm) 

- 2yn{tm)] + ^tCi{n,t^), (22) 

- 2x„(t„0] + Aty„(t„), (23) 

where tm — mAt. For the boundary points n = I and 
n = we use free boundary conditions, i.e., we hold 
xq = xi, uo = yi, XN = XN+i and ypf = yN+i for all 
times tm- We choose At — 0.1 in our simulations [23], 
and used chain lengths of size N = 1000 - 10000. The 
variable Ci('^j^m) is an independent Gaussian variable 
for all n and tm and distributed with zero mean and unit 
variance. 

The persistence probability for x„ upto time tm was 
obtained by keeping track of the fraction of a;„'s that have 
sgn[xn{tm)] same as sgn[xn{l)] for all times starting from 
1 to tm- The data are shown in fig. [TJ Typically each 
data curve in fig. [1] was obtained by averaging over 500 
thermal histories. We find that Pa{t) decays as a power 
law ~ t-^ with 9 ~ 0.360 ± 0.001. The latter value is in 
good agreement with the IIA estimate in Eq. (HJ). 



VI. SIMULATION OF THE GAUSSIAN 
PROCESS 

Since a Gaussian stationary process X{T) is com- 
pletely specified by its stationary correlator A{T), 
one can simulate the process by constructing a time- 
series with the same correlator. In the frequency do- 
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FIG. 2: The graph (dashed hne) shows P(T) vs. T from the 
simulation of the Gaussian process X{T). The upper short 
thick hne goes as 0.7 exp(~0.355r). 

main (Fourier space) the corresponding correlator is 

{X{uji)X{uj2)) = 27rJL(a;i)(5(u;i + W2), where A{uj) = 
/ dte''^'^A{T) is the Fourier transform of A{T). The lat- 
ter formula allows us to easily generate stochastic pro- 
cesses 

X{Lu)=f^{Lu)^Aso{io), (24) 

where fj{uj) is a Gaussian white noise with {fj{uji)fj(uj2)) — 
2tt6{uji + W2)). 

We performed simulations following the above route by 
first constructing random functions X{uj) as per Eq. ([M]) 
for discrete w's. Then we did a discrete inverse Fourier 
transform to obtain the times series X{T) 27]. After 
generating 10^ such random time-series oi X{T), we used 
them to calculate the probability density function P{T) 
of intervals between two consecutive zero-crossings. In 
the calculation, the time-step size used was ST — 1, and 



Tmax = 50 as we found that A{T) almost vanishes for 
T > 50. As stated earlier, in terms of the variable T, 
both P{T) and Po(r) decay as ~ exp(-6'T). Hence, 
from the decay of P{T) ~ exp(— ^T), we estimated 9. In 
Fig. [21 we have shown P{T) versus T, and we find the 
decay constant 6 w 0.355. The latter value is slightly 
smaller than the obtained from IIA and the Langevin 
simulation, because the step-size ST = 1 was a bit large 
and we missed some intervals smaller than that. 



VII. CONCLUSION 

In summary, we have studied the persistence probabil- 
ity of the X coordinate of a bead in the bulk of a Rouse 
polymer chain advected by a shear flow field. We have 
shown that the persistence probability decays as a power 
law in time at late times and the associated persistence 
exponent ^ « 0.36 is nontrivial. We have computed this 
exponent analytically within an independent interval ap- 
proximation and also determined it numerically by two 
different methods. The analytical result is in excellent 
agreement with the numerical simulations. 

There are several directions in which our work can be 
extended. Here we have considered the Rouse chain em- 
bedded in two spatial dimensions. It should be relatively 
straightforward to extend our method to calculate the 
persistence properties of the Rouse chain in higher di- 
mensions in presence of a transverse shear flow. It would 
also be of interest to study the persistence properties of 
the polymer chain in a more realistic setting going beyond 
the simple Rouse model, e.g, in presence of excluded vol- 
ume interactions. 
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